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NUMERICAL  METHODS  IN  GASDYNAMICS 


A.  A.  Dorodnitsyn 

1.  Introduction 

The  advent  of  high-speed  computing  machines  is  of  particularly 
important  value  in  numerical  methods  in  all  fields  of  the  exact  sciences 
among  them,  gasdynamics,  since  it  allows  us  to  obtain  solutions  of 
complete,  unsimplified  equations  with  an  accuracy  that  comple tely 
satisfies  and  even  exceeds  practical  requirements.  The  results  have, 
in  many  cases,  provided  much  more  exact  data  than  a  simulated  experiment 
to  say  nothing  of  the  fact  that  the  calculation  requires  much  lower 
expenditures  than  does  an  experiment. 

The  characteristic  feature  of  gasdynamics  processes  is  that, 
generally  speaking,  they  are  accompanied  by  discontinuities  in  the 
functions  defining  the  process  (pressure,  density,  velocity,  etc.). 
Therefore,  for  the  numerical  solution  of  problems  in  gasdynamics  it 
was  necessary  to  develop  mathematical  methods  which  allowed  us  to 
obtain  discontinuous  solutions. 

It  often  happens  that  a  method  of  solution  which  is  strictly 
valued  for  a  certain  class  of  functions  is  actually  applicable  to  a 
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a  considerably  broader  class.  Comparatively  recently,  the  hope  still 
existed  that  the  methods  of  solution  of  partial  differential  equations 
developed  for  continuous  and  smooth  functions  would  prove  effective  for 
discontinuous  solutions  or  equation  coefficients.  However,  these 
hopes  have  not  been  justified.  I  shall  give  one  very  simple  example 
( suggested  by  A.  A.  Samarskiy) ,  which  will  illustrate  how  a  method 
which  is  very  good  for  continuous  and  smooth  functions  leads  to 
completely  incorrect  results  when  it  is  applied  to  discontinuous 
functions. 

Let  us  consider  the  simplest  one-dimensional  equation 

(1,1) 

(it  may  be  treated,  for  example,  as  a  heat-conductivity  equation  or 
a  one-dimensional  discontinuity  equation)  with  boundary  conditions 
(1.2)  “(°) 0,  u(l) »  i.. 

For  continuous  and  smooth  k(x)  very  good  numerical  solutions  are 
given  by  the  finite -difference  scheme 

(1.5)  + 


(h  is  the  network  spacing  along  x) . 

Now  let  k(x)  be  a  discontinuous  function 


(1.4) 


k{  —  const 


where 

where 


An  accurate  solution  of  (1.1)  satisfying  the  condition  of  continu¬ 
ity  of  u(x)  ( tenjperature)  and  of  flow  k  du/dx  =  0  (let  us  say  of  heat 
flow)  is  easily  found 
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(1.5) 


’  ^ 

T+(^l)f  when  *<*• 
(x-i)f+j?  when  x>t 


An  exact  solution  of  finite-difference  equation  (1.4)  is  also 
found  without  difficulty.  In  fact,  let  the  point  of  discontinuity 
4  lie  between  x^  and  Xj  +  ^  [if  the  entire  interval  (0,1)  is  broken 
up  into  n  spacings  h  =  1/n,  j/n  <  4  <  ( J  +  l)/n],  when  i  ^  J 

«,  =  Cx, 

(the  constant  C  will  be  determined  later  on  from  the  condition  * 
Then,  checking  the  two  values  Uj  +  ani  u^  +  2,  we  obtain 

L  (x+3)f5x-l)*J 


When  i  >  J  +  2  Eq.  (1.5)  again  reduces  to 

«hi-2«,+“,-i  =  0 

and  when  i  >  J  +  2  elementary  calculation  gives 

M  r[(S-")(3*  +  ') v,/,  (5-,)(3x+_l)\  ,  2(x-l)(13x+7)  1. 

(1.7)  M'  rUvf.*\5.v  ■n1,T  (x+3)(5x-l)  "J 

Setting  i  =  n,  and  taking  into  account  the  fact  that  u^ 
we  obtain  for  C  the  expression; 

Q  _ _ _ _ _ _ 1 _ . 

|~ (5  **)  (3*+l)  /.  (5-x)(3x+l)\  2(x—l)(13w+7)  L]  * 

L(x+3)(5x-1)+(  (x+3)(5x-l)/*'+l4'  (x+3)f5*-l;  *J 

By  passage  to  the  limit  h -»  0  (x^  +  ^  ->  4)  we  obtain: 


C*+3)  (5y  n'1  *> 
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and  for  u 


(1.8) 


X 


M- 


(5-  x)(3x+l) . 
+3)  (5x— l)'1 


£  + 
S  + 


(5 -x)  (3x4-1) 
(x 4- 3)  (5x—  1) 
(5-x)  (3x4-1) 
(x-f-3)  (5x— lj 


(*~t) 

I 

0-0 


wheri*^** 


when  x>f. 


Eqs.  (1.5)  and  (1.8)  coincide  only  when  x  «  1  (that  is,  ki  =  k2)  , 
thus  the  solution  of  finite-difference  equation  (1.5),  even  though, 
generally  speaking,  it  approaches  a  definite  limit  (1.8)  [only  for 
certain  values  of  £  and  x  can  the  denominator  in  Eqs.  (1.8)  become 
zero],  does  not  give  the  solution  to  the  stated  problem.  As  is  easily 
seen  there  corresponds  to  the  solution  of  the  finite-difference  scheme 
the  presence  of  a  heat  source  at  the  point  x  «  £ 


3A,(x-l)» 

(5—f)(3*+ 1)4  8(**—  l)f  I 


i.e.,  the  finite -difference  scheme  does  not  ensure  the  integral  condi¬ 
tion  of  conservation  of  heat  flow  at  points  of  discontinuity. 

It  should  be  noted  that  the  property  of  finite-difference  schemes 
observed  in  this  example  forces  us  to  pay  very  close  attention  to  their 
use.  In  the  solution  of  nonlinear  gasdynamics  equations,  we  are  seldom 
successful  in  rigidly  proving  that  the  method  converges  to  the  required 
solution.  The  convergence  is  usually  verified  "empirically, "  by 
carrying  out  a  series  of  calculations  with  ever  decreasing  network 
spacings,  and  if,  beginning  with  some  value  of  h  a  further  decrease 
in  it  do«s  not  change  the  result  within  the  prescribed  accuracy  of 
the  calculation,  we  usually  accept  the  convergence  "on  faith". 

The  example  Just  presented  shows  that  in  the  case  of  discontinuous 


solutions  this  criterion  of  the  correctness  of  the  method  leads  to 


error. 

The  source  of  error  is  the  fact  that  the  finite  difference  scheme 
chosen  for  the  solution  of  the  equation  does  not  approximate  the  inte.- 
gral  law  of  conservation  of  heat  flow  in  the  presence  of  discontinuities. 
Thus  we  must  construct  approximating  operators  so  that  in  the  limiting 
case  of  discontinuity  they  will  accurately  represent  the  integral  laws 
of  conservation. 

2,  Method  of  Finite  Differences 

The  method  of  finite  difference  is  the  most  common  for  the  solution 
of  partial  differential  equations.  First  developed  for  equations  of 
the  elliptic  type,  it  was  then  generalized  for  hyperbolic  and  parabolic 
equations.  In  recent  years  this  method  has  found  wide  application  to 
the  solution  of  nonlinear  partial  differential  equations,  as  well  as 
in  cases  when  their  solutions  are  discontinuous.  I  should  like  to 
dwell  in  detail  on  the  latter  problem. 

Two  methods  of  approaching  the  solution  of  this  problem  are  recog¬ 
nized  at  the  present  time.  The  first  method  dates  back  to  Richtmayer 
and  consists  of  the  introduction  of  an  "artificial  viscosity"  into 
the  Initial  equations.  Here,  therefore,  the  original  problem  is  not 
solved,  but  a  modified  physical  problem,  in  which  the  discontinuities 
are  excluded.  The  discontinuities  are  replaced  by  regions  of  abrupt 
change  in  values.  The  solution  of  the  original  problem  is  obtained  as 
the  limit  of  the  solution  to  the  modified  problem  as  the  viscosity 
coefficient  tends  to  zero.  The  second  method,  that  about  which  we 
were  speaking  in  the  Introduction,  consists  of  the  special  construction 
of  finite -difference  schemes  for  which  the  Integral  conservation  laws 
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are  satisfied  In  the  limiting  case  of  a  discontinuity. 

I  shall  dwell  only  briefly  on  the  artificial-viscosity  method, 
pointing  out  those  difficulties  which  are  encountered  in  its  practical 
application. 

In  the  numerical  solution  of  equations  by  the  artififial -viscosity 
method,  it  is  actually  necessary  to  pass  to  a  limit  twice:  first  for 
a  fixed  viscosity  coefficient  the  network  spacing  tends  to  zero  and 
then  (after  obtaining  a  number  of  solutions  with  different  viscosity 
coefficients)  the  viscosity  coefficient  tends  to  zero.  Of  course, 
this  method  of  calculation  will  be  extremely  cumbersome.  In  reality 
the  calculation  is  carried  out  in  such  a  way  that  the  coefficient  of 
viscosity  decreases  simultaneously  with  the  network  spacing  ( in  practice 
the  viscosity  coefficient  is  chosen  so  that  the  shock  wave  will  be 
"spread"  over  5-6  network  spacings)  . 

The  introduction  of  viscosity  generally  lowers  the  accuracy  of  the 
calculation.  However,  the  viscosity  is  needed  only  where  there  are 
shock  waves  and  is  not  needed  where  the  solution  is  smooth.  Therefore, 
the  viscosity  coefficient  will  not  be  taken  constant,  but  as  a  function 
of  the  velocity  gradient,  in  the  following  form. 

If  the  initial  equation  was 

(2.i) 

[u  and  F(u)  may  be  considered  vectors],  the  modified  equation  may  be 
written : 


(2.2) 


dtL  4.  dFM  _ ,  A I  w(d±  \  \ 

dt  **"•  bx  <)*[  \  dx/ dJC I’ 


where  is  a  positive  function  which  increases  with  an  increase 

in  the  velocity  gradient  du/dx  [we  may  also  set  1(0)  *  0].  Such  a 
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device  accomplishes  a  decrease  in  the  influence  of  viscosity  where 
its  introduction  is  disadvantageous. 

The  other  difficulty  is  the  fact  that  in  replacing  Eq.  (2.2)  by 
a  finite-difference  equation,  convergence  is  ensured  in  the  presence 
of  discontinuities  and  when  "e-*  0  only  when  this  difference  equation 
has  first-order  accuracy.  Therefore,  in  order  to  obtain  a  sufficiently 
accurate  result,  it  is  necessary  to  carry  out  the  calculation  with  a 
very  large  number  of  network  nodes.  The  calculation  becomes  exceedingly 
cumbersome  even  for  the  best  high-speed  computing  machines. 

Of  course  convergence  of  difference  schemes  of  a  higher  order 
of  accuracy  is  disrupted  in  the  vicinity  of  discontinuities.  In  regions 
of  smooth  solution,  finite-difference  schemes  of  a  higher  order  of 
accuracy  improve  the  convergence.  In  order  to  ensure  convergence  and, 
at  the  same  time,  not  impair  the  convergnece  outside  the  vicinity  of 
the  shock  wave,  finite -difference  equations  are  used,  which  may  condi¬ 
tionally  be  called  "schemes  of  intermediate  accuracy".  Let  x)(I>(w)~o 
be  a  finite-difference  equation  of  first-order  accuracy,  approximating 
Eq.  (2.2),  and  .(**'(")  =  0-  a  finite  difference  equation  of  second- 
order  accuracy.  Any  equation  of  the  form 

(2.3)  «(?"'(«)+(!  —  =  () 

will  also  approximate  Eq.  (2.2)  and  if  a  >  0  and  fixed,  the  expression 
on  the  right  side  of  (2.3)  is  of  first-order  accuracy.  But  it  is 
possible  to  take  the  quantity  a  as  a  variable;  to  take  it  as  equal 
to  or  close  to  unity  where  there  is  "danger”  of  the  formation  of  shock 
waves,  and,  conversely,  where  the  solution  is  smooth  to  take  it  as 
small  or  equal  to  zero.  Of  course,  this  process  of  selecting  a  must 
be  automated  for  machine  calculation,  for  which  it  is  necessary  to 
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compute  the  "criterion  of  wave  generation"  by  the  values  of  u  at  a 
number  of  neighboring  points  and  to  assign  a  as  some  function  of  this 
criterion.  This  method  allows  us  to  increase  the  accuracy  of  the 
calculation.  However,  it  leads  to  a  very  appreciable  complication  of 
the  logical  scheme  of  the  calculation,  that  is,  it  significantly  compli¬ 
cates  the  programming  of  problems  for  high-speed  computing  machines. 

Let  us  now  consider  the  second  method  of  solving  discontinuous 
problems.  The  differential  equations  of  gasdynamics  are  expressions 
of  the  laws  of  the  conservation  of  mass,  momentum,  and  energy.  There¬ 
fore  we  may  list  them  in  "divergent"  form. 

In  the  case  of  two  independent  variables  (two  coordinates  or 
one  spacial  variable  and  time) ,  the  equations  may  be  written  in 
"divergent"  form  as  follows: 


(2.4) 


where  P,  Q,  and  F  are  some  functions  of  unknown  quantities  and,  perhaps, 
independent  variables.  Right  sides  which  differ  from  zero  may  occur 
in  the  presence  of  sources  (for  example  heat  liberation  due  to  chemical 
reactions).  Integrating  Eq.  (2.4)  over  some  domain  bounded  by  the 
contour  2:  ,  we  arrive  at  the  integral  relation: 

(2.5)  J/Vv  -o,/.v  fjFdxdy, 

•'  '■/> 

which  will  be  valid  not  only  for  continuous  P  and  Q,  but  for  piecewise- 
continuous  as  well.  Thus  equations  in  the  form  of  (2.5)  will  always 
be  effective  for  gasdynamic  processes,  while  Eqs.  (2.4)  may  lose  meaning 
in  the  vicinity  of  discontinuities.  Specifically,  all  the  conditions 
for  the  discontinuities  are  obtained  from  Eq.  (2.5) . 

Therefore,  even  when  constructing  numerical  methods  of  solution 
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suitable  for  application  to  the  case  of  discontinuities,  it  is  necessary 
to  proceed  from  equations  in  the  form  of  (2.5),  and  not  from  partial 
differential  equations. 

For  clarity  we  shall  examine  the  construction  of  finite -difference 
equations  for  one -dimensional  plane  motion  of  a  gas.  By  using  the 
Lagrangian  coordinate  system  in  which  the  variables  are  the  time  Jb  and 
the  mass  of  a  gas  column  of  unit  area  and  cross  section  x  (read  from 
some  initial  point) .  In  this  coordinate  system,  the  equations  of  the 


problem  are  written  in  the  form: 


(2.6) 


dt  +  dx  ’  dt 


(here  u  is  the  velocity;  £  the  pressure;  v  the  specific  volume;  aid 
E  the  internal  energy,  for  an  ideal  gas  £  =-/>r/(x-l),  ,  where  x  =  Cp/C. 

is  the  ratio  of  the  specific  heats). 

Corresponding  to  these  differential  equations  are  the  integral 
relations: 

(2.7)  / \udx-pdt)  —  0,  j  ( vdx+udt )  =  0,  /[(£+t)*-H=0- 

After  taking  a  network  rectangle  with  vertices 

(*»i  f«)i  (*mtl  >  0>  (*»+l>  f*+0>  (*«.• 

over  the  contour.^,  the  integral  relations  (2.7)  may  be  presented  in 
the  form: 


(2.8) 


--  (i. 


Here  // -=xm,,  T  ,  the  subscript  m  denotes  the  value  of 

the  function  when  x  =  xm,  the  subscript  n  the  value  of  the  function 
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whan  t  »  tft,  the  subscript  m  +  1/2  some  mean  value  in  the  interval 
(xm,  xm  i )  and  n  +  1/2  some  mean  value  in  the  interval  (tn,  tn  +  j.) . 

If  all  quantities  were  continuous,  we  could  easily  relate  the  mean 
values  of  the  quantities  to  the  values  at  the  nodes  of  the  network 
(for  example,  by  taking  the  half-sum  of  neighboring  nodal  values). 
However,  since  we  wish  to  construct  a  method  of  calculation  which  is 
also  suitable  in  the  presence  of  shock  waves,  we  should  consider  any 
point  in  the  region  a  possible  point  of  discontinuity.  In  the  finite- 
difference  interpretation,  this  means  every  network  node  should  be 
considered  a  location  of  discontinuity.  Thus  at  the  moment  of  time  tR 
we  consider  the  point  xm  a  point  of  discontinuity  to  the  right  of  which 
the  values  of  S  are  equal  to  +  ,  and  to  the  left,  equal  to 

Sm  -  1/2  *  Then  waves  ( shock  or  rarification)  will  leave  in  both 
directions  from  thepoint  xffl,  at  that  same  point  xm  values  of  the  quan¬ 
tities  will  be  established  corresponding  to  the  integral  conditions  for 
the  discontinuity  [*«  -  p\  =  o,  [an  «]  =  0,  [c(£+  u'/l)  ~pu]  -  0,  where  c  is  the 

velocity  of  the  shock  wave.  This  value  will  be  retained  until  waves 
from  other  network  nodes  reach  this  point.  But  if  the  time  spacing  t 
is  chosen  so  that  this  does  not  occur  (this  condition  also  satisfies 
the  convergence  condition  of  the  finite-difference  method  for  hyperbolic 
equations),  the  quantities  Pm,  vm  etc.  will  remain  constant  throughout 
the  entire  interval. 

The  dependences  between  Pm,  Pm  +  1/2*  and  Pm  -  1/2  may  be  Presented 
(after  transforming  the  conditions  for  discontinuities  In  the  form 

(2.9)  f(:“  2.'  i 

.  pJ,c*  +  cn)  =  c„Pm  1  +  rA  Pm+ 1  +  cA  cH(um 
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Here  c,  and  c..  are  the  absolute  values  of  the  velocities  of  the  left 
l  r 

and  rigit  shock  waves  respectively.  Eqs.  (2.9)  are  valid  if  both  shock 
waves  p*>pm  ft-  .  For  rarification  waves,  Eqs.  (2.9)  should 

be  replaced  by  the  corresponding  conditions  for  rarification  waves. 

After  determining  Pm  from  Eq.  (2.9),  vm  is  found  elementarily 

f'n.(c«  +  c*)  ~  ClU* i-J  +^“"1  —  Pm+\ 

and  Eqs.  (2.8)  allow  us  to  make  the  time  spacing. 

The  method  set  forth  above  was  proposed  by  S.  K.  Godunov  in  1952 
and  has  been  successfully  applied  to  the  solution  of  practical  problems. 

In  constructing  a  finite-difference  approximation  for  partial 
differential  equations  to  obtaining  discontinuous  solutions,  each 
network  node  must  be  considered  a  possible  point  of  discontinuity.  The 
relationship  between  the  quantities  for  the  discontinuity  must  be 
taken  in  accordance  with  the  integral  conditions  of  conservation  (i.e. 
to  take  these  relations  as  exact,  as  in  the  example  set  forth  above, 
or  approximate,  but  such  as  to  ensure  the  prescribed  accuracy  of 
calculation) . 


3.  Method  of  Integral  Relations 

A  very  effective  method  for  solving  problems  on  high-speed  computing 
machines  is  the  method  of  approximate  reduction  of  partial  differential 
equations  to  systems  of  ordinary  differential  equations.  This  method 
allows  us  to  make  use  of  the  highly  developed  apparatus  for  the  solution 
of  ordinary  differential  equations,  and,  in  addition,  in  the  presence 
of  unbounded  domains  it  allows  us  to  use  other  well  developed  methods 
for  the  asymptotic  solution  of  ordinary  differential  equations. 

Difficulties  arise  in  the  application  of  this  method  (especially 
with  nonlinear  equations)  when  it  is  necessary  to  solve  a  boundary- 
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value  problem  for  the  approximate  systems  of  ordinary  differential 
equations.  Therefore,  this  method  is  effective  for  boundary-value 
problems  in  those  cases  which  allow  us  to  find  a  sufficiently  accurate 
solution  when  the  order  of  the  approximate  system  of  ordinary  equations 
is  low.  As  experience  has  shown,  a  highly  satisfactory  accuracy  is 
attained  when  integral  methods  are  used,  i.e..  Integral  methods  when 
the  integral  relations  expressing  the  laws  of  conservation  are  set  as 
the  basis  for  construction  of  an  approximate  system  of  ordinary  differ¬ 
ential  equations. 

As  stated  above,  the  gasdynamics  equations  may  be  presented  in 
"divergent"  form: 

^  dxp^x,y,Ul . M")+  ly^x,y> •••*“") ” y' M* . “») 

(*  =  1.2, 


where  ux,  ...,  are  unknown  functions. 

Let  it  be  necessary  to  solve  system  of  equations  (3.1)  in  domain 
,0 ,  which  has  the  shape  of  a  curvilinear  rectangle,  the  boundaries 
of  which,  let  us  say,  are  the  lines 
(3.2)  v =- a,  .v  =  b,  y=o,^  =  a(*) 

(in  special  cases  a,  may  be  equal  to  -  •»,  b  *  +») .  Dividing  the  domain 
up  into  N  strips  of  width  6/H  and  integrating  each  of  the  equations  of 
system  (3.1)  with  respect  to  £  across  each  strip,  we  obtain  the  system 
of  integral  relations 


(5.5) 


(  »*+i 

~dx  J  *+i~  Qi,k  —  j' Ftdy, 

•*  »* 

*—  1,  2 . tv,  A  =  0,  1 . N—  1;  yk=  ^i. 
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The  upper  boundary  of  the  domain  y  -  d(x)  can  be  a  given  function, 
and  it  can  also  be  an  unknown  function,  subject  to  determination.  In 
the  first  case  a  total  of  n  boundary  conditions  should  be  given  at  the 
boundaries  y  *  0  and  y  «  6(x),  and  in  the  second  case,  n  +  1  boundary 
conditions  should  be  given. 

Of  course  the  lower  boundary  of  the  domain  may  also  be  curvilinear. 
The  generalization  of  integral  relations  (5.3)  for  this  case  is  obvious. 

Integral  relations  (3. 3)  present  themselves  physically  as  the 
conservation  laws  (mass,  energy,  momentum,  etc.)  set  down  for  the  strips. 

If  we  now  represent  the  functions  and  with  the  aid  of  some 

interpolation  formulas  by  their  values  k  and  F^  k  at  the  boundaries 

of  the  strips,  the  integrals  in  the  relations  (3.3)  are  approximately 
represented  in  the  form 


(3.4) 


/■ 


»*+i 

P,dy  =  d  yAk  tp L, 


»— 0 


and  analogously  for  the  integral  of  P1<  In  Eq.  (3.4)  the  coefficients 

A  are  definite  numbers  which  depend  on  the  chosen  interpolation 
k,  v 

formula. 

The  substitution  of  Eq.  (3.4)  into  integral  relations  (3.3)  leads 
to  a  system  of  nN  ordinary  differential  equations  relative  to  n(N  +  1) 
unknown  functions  u^  ^  (or  n(N  +  1)  +  1  unknown  functions  if  6(x)  is 
also  previously  unknown) .  The  addition  to  the  system  of  n  (or  n  +  1 
when  6(x)  is  not  given  )  conditions  for  the  boundaries  y  =  0  and  y  =  6(x) 
closes  the  system,  and  thus  the  problem  is  reduced  to  a  closed  system 
of  ordinary  differential  equations. 

Note  that  the  selection  of  the  system  of  functions,  with  which 
the  interpolation  formulas  for  P.^  and  F^^  are  constructed,  has  signifi¬ 
cant  importance  for  the  accuracy  of  the  calculation.  The  "inaptness" 
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of  the  selection  of  the  system  of  interpolation  functions  depends,  of 
course,  on  how  well  we  represent  qualitative  behavior  of  the  solution. 

As  is  apparent  from  the  above, 

characteristic  the  method  of  integral  relations  is 

easily  applied  to  the  case  of  unknown 

boundaries  of  a  domain.  This  has 

special  significance  in  gasdynamlcs 

problems,  where  the  shape  of  the 

shock  wave  and  the  boundaries  of 

Fig.  1.  Flow-past  diagram  the  domains  of  influence  are  not 
at  \  *  1. 

known  beforehand. 

Let  us  consider,  for  example,  the  problem  of  flow  past  a  body 
with  a  velocity  exactly  equal  to  the  speed  of  sound.  This  problem 
was  solved  by  P.  I.  Chushkin.  For  simplicity,  we  3hall  consider  a 
symmetrical  body.  Here  it  is  necessary  to  find  a  solution  to  the 
equations  of  gasdynamlcs  in  the  domain  bounded  by  the  axis  of  symmetry, 
the  body  contour,  and  the  characteristic  of  the  first  family  tangent 
to  the  sound  curve  at  infinity  (Fig.  1) .  The  shape  of  this  character¬ 
istic  is  unknown.  Corresponding  to  the  line  y  =  0  in  the  general 
theory  of  the  method  is  the  axis  of  symmetry;  corresponding  to  the 
curve  y  =  6(x)  is  the  above-mentioned  characteristic.  We  have  two 
z'elations  for  this  characteristic:  the  differential  equation  of  the 
characteristic  itself  and  the  relation  between  the  increments  in  the 
dope  of  the  velocity  vector  and  the  increments  in  the  velocity  itself. 

In  this  problem  the  number  of  unknown  functions  n  is  equal  to  2 
(let  us  say  the  velocity  v  and  the  slope  of  the  velocity  *) .  For  the 
axis  of  symmetry  we  have  one  condition  (*  =  0) ,  and  2  conditions  for 
the  unknown  characteristic,  a  total  therefore  of  3  =  n  +  1  conditions. 
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In  the  problem  of  the  supersonic  flow  past  a  blunt  body,  i.e., 
with  a  detached  shock  wave,  the  shape  of  the  shock  wave  is  unknown. 
Corresponding  to  the  shock  wave  Is  the  boundary  y  -  6(x)  of  the  general 
theory;  corresponding  to  the  body  contour  is  the  boundary  y  *»  0  (Pig.  2) . 
In  this  problem  there  are  a  total  of  4  unknown  functions  (e.g.,  the 
two  components  of  the  velocity  u  and  v,  the  density  • ,  and  the  pressure 
£  ) .  The  condition  of  no  flow  by  the  body  contour  provides  one  boundary 
condition,  and  for  the  shock  wave  we  have  four  conditions  (each  of  the 
quantities  u,  v,  *,  £  expressed  in  terms  of  the  slope  of  the  shock 
wave).  Altogether,  therefore,  we  have  5  ■  n  +  1  conditions,  i.e., 
as  many  as  are  necessary  for  the  construction  of  a  closed  system  of 
ordinary  differential  equations,  in  which  the  width  of  the  region  will 
enter  as  one  of  the  unknown  functions.  The  problem  of  the  supersonic 
flow  past  bodies  with  a  detached  shock  wave  was  solved  by  0.  M. 
Belotserkovskiy. 

shock  wave  The  two  P1*0*516®8  given  here, 

as  far  as  I  know,  have  not  been 

solved  by  any  other  method. 

It  is  Interesting  to  Illustrate 

the  actual  rapidity  of  convergence 

of  the  method  of  Integral  relations 

Pig.  2.  Flow-past  diagram  using  the  results  of  P.  I.  Chuskin 
with  detached  shock  wave. 

and  0.  M.  Belotserkovskiy.  Calcula¬ 
tions  of  the  critical  Mach  number  for  ellipses  and  ellipsoids  according 
In  the  first,  second,  and  third  approximations  are  presented  in  Pig. 

5.  The  velocity  distribution  around  a  circle  when  M  =  M„_  in  the 
presence  of  circulation  is  given  in  Fig.  4. 
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fig.  5.  Critical  Mach  a unbar 
distribution  for  ellipses  and 
ellipsoids. 


fig.  3,  Pressure  distribution 
around  a  circla  whan  1^*  L 


fig.  4*  Velocity  distribution 
around  a  oircla  with  circula¬ 
tion  at  M  . 

cr 


istie  and  along  tbs  axis  for 
a  circla  whan  hM  *  1. 


around  a  circle  whan  M  *  5. 

fig.  7.  Shape  of  tba  shock 
ways  for  a  circla  that  M  *  J. 

w 
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In  Fig.  5  and  6  are  listed  the  results  of  calculations  of  the 
flow  past  a  circle  of  a  gas  stream  at  the  speed  of  sound  at  infinity. 

The  shape  of  the  shock  wave  is  shown  in  Fig.  7,  and  the  pressure 
distribution  around  a  cylinder  with  flow  past  by  a  supersonic  stream 
is  shown  in  Fig.  8. 

As  is  apparent  from  the  graphs,  the  calculation  is  sufficiently 
accurate  in  the  second  approximation,  in  which  the  entire  domain  has 
been  subdivided  into  only  two  strips,  we  can  hardly  say  here  that  the 
law  of  the  convergence  to  zero  of  the  error  of  the  approximation 
appeared  in  the  second  or  third  approximation.  From  a  practical  point 
of  view,  the  limiting  law  of  convergence  to  zero  of  the  error  of  the 
method  is  of  little  interest.  It  is  important  for  the  method  that  the 
accuracy  needed  in  practice  be  attained  in  low-order  approximations. 

And  the  examples  presented  here  illustrate  the  effectiveness  of  the 
method  of  integral  relations  in  this  respect. 

4.  Remark  on  the  Method  of  Characteristics 

The  classical  problems  of  the  method  of  characteristics  •— 

Cauchy's  problem,  Ooursat's  problem,  and  also  the  determination  of  the 
gas  flow  along  a  given  characteristic  curve  of  the  first  family  and 
the  contour  of  the  body  past  which  the  flow  occurs  —  at  the  present 
time  (with  calculation  on  high-speed  computing  machines)  can  be 
considered  extremely  simple  problems. 

Cases  in  which  characteristic  curves  of  one  family  intersect,  i.e. 
in  which  shock  waves  arise,  present  the  greatest  difficulty  when  the 
method  of  characteristics  is  used. 

In  steady-state  problems  of  aerodynamics,  it  is  possible  in  most 
cases  to  determine  the  presence  of  shock  waves  and  their  approximate 
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position  beforehand.  This  is  considerably  harder  to  do  in  nonsteady- 
state  problems.  Therefore  the  method  of  finite  differences  is  chiefly 
used  in  the  solution  of  nonsteady-state  problems,  though  the  convergence 
of  the  method  of  characteristics  is  considerably  better  than  that  of 
the  method  of  finite  differences. 

Known  difficulties  are  encountered  when  the  method  of  character¬ 
istics  is  used  in  the  vicinity  of  the  sound  curve  and  when  constructing 
a  head-attached  shock  wave. 

In  the  first  case  it  is  most  advisable  to  depart  from  the  sound 
curve  with  the  aid  of  series  (it  is,  in  practice,  sufficient  to  use  one 
or  two  terms  of  the  series) . 

In  calculating  by  the  method  of  characteristics  with  simultaneous 
construction  of  the  shock  wave,  it  is  difficult  to  ensure  stability  of 
the  calculatl on,  especially  if  the  angles  of  intersection  of  the 
characteristic  curves  of  the  first  family  with  the  shock  wave  are 
small.  Here  it  is  also  more  advisable  to  use  series  in  the  construction 
of  the  initial  element  of  the  shock  wave. 

5.  Conclusion 

In  conclusion,  I  wish  to  note  that  from  my  point  of  view  the  most 
important  problem  at  the  present  time  in  the  field  of  the  numerical 
methods  in  gasdynamics  is  the  development  of  effective  methods  for  the 
solution  of  three-dimensional  problems  (spatial  steady-state  flows  and 
nonsteady-state  two-dimensional  ones) . 

Of  course,  one  cam  not  say  that  all  two-dimensional  problems  have 
already  been  solved  and  that  nothing  remains  to  be  done  with  regard  to 
improving  the  methods  of  calculation.  However,  one  can  say  that  with 
high-speed  computing  machines  almost  any  two-dimensional  problem  can 
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be  solved,  even  If  in  a  very  cumbersome  way.  Although  only  relatively 
simple  three-dimensional  problems  are  being  solved,  there  still  remains 
in  this  problem  many  fundamental  difficulties,  and  not  only  confutation 
difficulties.  As  regards  the  practical  significance  of  three-dimen¬ 
sional  problems,  this  is  obvious  to  everyone. 

The  development  of  methods  for  the  numerical  solution  of  three- 
dimensional  problems  would  give  scientists  and  engineers  in  the  field 
of  aerodynamics  a  powerful  tool  in  their  scientific  research,  as  well 
as  in  the  solution  of  actual  engineering  problems. 
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